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1.  Introduction 


The  continuing  development  of  microelectromechanical  systems  is  pointing  to  the  possibility  of 
mounting  complete  sensor  systems  on  medium-  and  small-caliber  projectiles  as  part  of  an  actively 
controlled  smart  munition.  Two  important  technical  challenges  in  achieving  this  goal  are  the 
development  of  small,  rugged  sensor  suites,  and  control  mechanisms.  With  regard  to  the 
development  of  control  mechanisms,  several  concepts  have  emerged  that  produce  controllable 
impulsive  lateral  forces  on  a  projectile  body.  For  example,  Harkins  and  Brown  (7)  considered  the 
use  of  a  set  of  lateral  pulse  jets  or  squibs  to  reduce  dispersion  of  a  rocket  by  firing  squibs  to  minimize 
projectile  pitch  and  yaw  rate.  For  the  notional  cases  evaluated,  dispersion  was  reduced  by  a  factor  of 
5.  Jitpraphai  and  Costello  (2)  considered  the  same  type  of  control  mechanism  and  used  a  trajectory 
tracking  flight  control  system  to  improve  impact  point  performance  of  a  direct  fire  rocket  equipped 
with  a  ring  of  squibs.  Amitay  et  al.  (3)  considered  the  use  of  synthetic  jet  actuators  as  a  control 
mechanism  on  lifting  bodies.  For  a  spinning  projectile  incorporating  a  synthetic  jet  actuator  for 
control,  the  synthetic  jet  actuator  is  activated  over  a  small  portion  of  a  roll  cycle  leading  to  a  train  of 
lateral  pulse  forces  acting  on  the  projectile  near  the  synthetic  jet  actuator  cavity. 

The  design  of  flight  control  systems  for  fin-stabilized  configurations  is  well  established  in  the  missile 
community  and  the  control  response  to  force  and  moment  inputs  is  reasonably  well  understood. 
Generally  these  configurations  are  treated  largely  in  the  same  manner  as  airplanes.  While  the 
uncontrolled  dynamics  of  spinning  projectiles,  both  fin  stabilized  and  spin  stabilized,  have  been 
extensively  studied  in  the  ballistics  community,  issues  with  regard  to  control  response  have  received 
considerably  less  attention  due  to  the  lack  of  practical  application  of  control  technology  to  spinning 
projectiles.  Using  projectile  linear  theory,  this  report  analytically  investigates  several  aspects  of  the 
response  of  a  spinning  projectile  to  lateral  pulse  forces  including  swerve  response  magnitude  and 
phase  angle,  impulse  force  smearing,  and  yaw  cycle  pulse  timing.  The  report  begins  with  a 
discussion  of  the  basic  projectile  dynamic  model  followed  by  judicious  simplifications  to  these 
equations  that  result  in  the  projectile  linear  theoiy  equations.  The  solution  of  the  projectile  linear 
theory  equations  are  used  to  shed  light  on  intuitive  and  subtle  factors  that  influence  swerve  response 
of  a  projectile  exposed  to  lateral  impulsive  loads. 


2.  Projectile  Dynamic  Model 


It  is  well  known  that  the  motion  of  most  projectile  configurations  can  be  captured  using  a  rigid 
body  6  degrees  of  freedom  dynamic  model  (4, 5).  The  degrees  of  freedom  include  three  position 
components  of  the  mass  center  of  the  projectile  as  well  as  three  Euler  orientation  angles  of  the 
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body.  Figures  1  and  2  show  schematics  of  the  dynamic  model  degrees  of  freedom.  The 
equations  of  motion  are  provided  in  equations  1—4. 


Figure  2.  Projectile  orientation  definition. 


u)  fX/ml  TO  -r  q  [u 
<  v  >  =  <Y /m  •-  r  0  -p  v»; 
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(1) 


(2) 


(3) 


The  applied  forces  in  the  body  frame  that  appear  in  equation  5  contain  contributions  from  weight 
(W),  air  loads  (A),  and  lateral  pulse  forces  (C). 


The  weight  force  resolved  into  projectile  body  coordinates  is  given  by  equation  6, 


(6) 


The  air  loads  are  split  in  two  components,  the  standard  aerodynamic  forces  and  the  Magnus 
forces, 

xA|  [x„l  [X™ 

■YA  =  Y  L+  Y  l  (7) 

J'A.  ,^AS.  .^AM. 
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In  equation  12,  Y,  and  Zr  represent  the  lateral  components  of  the  impulse.  The  total  magnitude 
of  the  impulse  on  the  projectile  is  Fj  =  ^/Y,2  +  Z\ .  The  function  d(t)  is  a  constant  positive  value 

when  the  pulse  is  active  and  is  zero  otherwise.  This  function  integrates  to  unity  for  a  single 
pulse. 

The  applied  moments  about  the  projectile  mass  center  are  due  to  aerodynamic  forces  and 
moments  (A)  as  well  as  pulse  forces  (C). 
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The  aerodynamic  moments  caused  by  standard  and  Magnus  air  loads  are  computed  with  a  cross 
product  between  the  distance  vector  from  the  mass  center  to  the  force  application  point  and  the 
force  itself.  An  unsteady  aerodynamic  damping  moment  is  also  present,  which  provides  a 
damping  source  for  angular  motion. 
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All  aerodynamic  coefficients  and  the  center  of  pressures  are  a  function  of  the  Mach  number  of 
the  projectile  mass  center.  The  dynamic  model  previously  described  is  highly  nonlinear  due  to 
both  three-dimensional  rotational  kinematics  expressions  and  the  presence  of  complex 
aerodynamic  forces.  The  applicability  of  the  equations  of  motion  previously  shown,  have  been 
validated  over  the  past  60  years  at  aeroballistic  ranges  throughout  the  world  (5). 


3.  Projectile  Linear  Theory 


Pressed  with  the  need  to  predict  the  trajectory  and  stability  of  a  ballistic  shell  so  that  useful 
performance  data  could  be  generated  with  primitive  computers,  early  ballisticians  vigorously 
investigated  mathematical  simplifications  to  the  equations  of  motion  of  a  projectile.  What 
emerged  over  time  was  a  set  of  simplified  and  solvable,  yet  accurate,  linear  differential 
equations,  which  today  is  commonly  termed  “projectile  linear  theory.” 


The  governing  equations  developed  previously  are  expressed  in  the  body  reference  frame.  In 
projectile  linear  theory,  the  lateral  translational  and  rotational  velocity  components  are 
transformed  to  a  nonrolling  reference  frame.  The  nonrolling  frame  or  so-called  fixed  plane 
frame  proceeds  with  only  precession  and  nutation  rotations  from  an  inertial  reference  frame. 
Components  of  linear  and  angular  body  velocities  in  the  fixed  plane  frame  can  be  computed  from 
the  body  frame  components  of  the  same  vector  through  a  single-axis  rotational  transformation. 
For  example,  the  body  frame  components  of  die  projectile  mass  center  velocity  are  transformed 
to  the  fixed  plane  frame  by 
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It  should  be  noted  that  the  ~  superscript  indicates  the  vector  components  are  described  in  the 
fixed  plane  reference  frame.  In  projectile  linear  theory,  a  change  of  variables  from  station  line 
velocity  component,  u ,  to  total  velocity,  V ,  is  performed.  Equations  16  and  17  relate  V  and  u 
and  their  derivatives. 

V  =  Vu2  +  v2  +  w2  =  -\/u2  +  v2  +  w2  ;  (16) 


du  dv  dw  du  _  dv  ~  dw 

dV  _  dt  dt  dt  _  dt  dt _ dt_ 

dt  V  V 


(17) 


A  change  of  variables  from  time,  t ,  to  dimensionless  arc  length,  s ,  is  also  made.  Equation  18, 
as  defined  by  Murphy  (4),  gives  the  dimensionless  arc  length, 


s  = 


(18) 


Equations  19  and  20  relate  time  and  arc  length  derivatives  of  a  given  quantity  £ .  Dotted  terms 
refer  to  time  derivatives,  and  primed  terms  denote  dimensionless  arc  length  derivatives. 
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In  projectile  linear  theory,  several  assumptions  regarding  the  relative  size  of  different  quantities 
are  made  to  simplify  the  analysis.  Euler  yaw  and  pitch  angles  are  small  so  that:  sin(0) «  6 , 
cos(0)  « 1 ,  sin(vj/)  «  v|/ ,  and  cos(v)/)  « 1  and  the  aerodynamic  angle  of  attack  is  small  so  that 
a  =  w/V  and  P  =  v/V .  The  projectile  is  mass  balanced  such  that  1^  =  Ixz  =  ^yz  =  0  and 
Izz  =  Iyy  =>  Iyy  =  ^zz  s  ly  •  The  projectile  is  aerodynamically  symmetric  such  that 
Cy0  =  Czo  =0.  Quantities  V  and  <j>  are  large  compared  to  0,\|/,q,r,v,  and  w ,  such  that 

products  of  small  quantities  and  their  derivatives  are  negligible.  Application  of  these 
assumptions  results  in  equations  21-29. 
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(30) 


These  equations  are  a  coupled  set  of  linear  differential  equations  except  for  the  fact  that  the  total 
velocity,  V ,  appears  in  the  coefficients  of  many  of  the  dynamic  equations. 


4.  Projectile  Linear  Theory  Solution 


Linear  theory  offers  physical  insight  into  the  projectile  flight  dynamics  because  closed  form 
solutions  can  be  readily  obtained  (4).  Using  the  assumption  that  V  changes  slowly  with  respect 
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to  the  other  variables,  it  is  thus  considered  to  be  constant,  V  «  V0,  when  it  appears  as  a 

coefficient  in  all  dynamic  equations  except  its  own.  Moreover,  pitch  attitude  of  the  projectile  is 
regarded  as  constant  in  the  velocity  equation,  uncoupling  the  velocity  equation  from  the  system. 
The  angle-of-attack  dynamics  or  epicyclic  motion  in  equation  29,  together  with  the  roll  dynamics 
in  equation  28  are  uncoupled  and  form  a  linear  system  of  differential  equations.  In  projectile 
linear  theory,  the  Magnus  force  in  equations  25  and  26  is  typically  regarded  as  small  in 
comparison  to  the  other  aerodynamic  forces  and  is  shown  only  for  completeness.  In  further 
manipulation  of  the  equations,  all  Magnus  forces  will  be  dropped.  Magnus  moments  will  be 
retained  however,  due  to  the  magnitude  amplification  resulting  from  the  cross  product  between 
Magnus  force  and  its  respective  moment  arm. 


The  solution  to  the  differential  equation  equation  27,  for  the  forward  velocity,  is 


(31) 


where 


av  = 


pSDCX0 

2m 


(32) 


and 


bv  =  gD0o .  (33) 

When  0O  =  0 ,  the  velocity  solution  reduces  to  the  familiar  exponential  decay  form  (4). 

The  roll  dynamic  equation  is  a  nonhomogeneous  linear  differential  equation  with  the  following 
solution: 
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In  order  to  develop  the  swerve  closed  form  solution,  the  epicyclic  equations  must  first  be  solved 
because  the  lateral  translation  and  rotational  velocity  components  are  contained  in  the  attitude 
differential  equations,  and  the  attitudes  are  contained  within  the  swerve  differential  equations. 

The  epicyclic  differential  equations  consist  of  a  set  of  four  coupled  nonhomogeneous  differential 
equations.  The  homogeneous  solution  is  easily  formed  using  the  free  vibration  modes  and  mode 
shapes,  and  the  results  are  given  by  the  intrinsic  complex  expression, 
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The  sum  and  product  of  sf  and  ss  give  the  relations  found  in  the  next  expression, 
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Two  more  simplifications  based  on  size  are  introduced.  First,  neglect  the  product  of  damping 
and  secondly,  the  product  AE  is  neglected  because  the  density  ratio  is  assumed  small.  A 
solution  may  now  be  obtained  for  both  the  fast  and  the  slow  damping  factors  and  turning  rates 
for  the  translational  and  rotational  velocities. 
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(42) 


<t>s 


5.  Pulse  Force  and  Moment  Conditions 


The  pulse  force  applied  to  the  projectile  is  taken  to  be  a  lateral  impulsive  force,  and  this  force  is 
due  to  an  actuator  attached  to  the  projectile  body  (equation  29).  For  this  investigation,  the  force 
actuator  is  modeled  as  a  scaled  square  wave  pulse  so  that  the  resulting  force  and  moment 
components  in  the  nonrolling  frame  are  as  follows: 
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for  a  square  wave  pulse  of  length  Ln  that  is  initiated  at  sn  represented  as 
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Note  that  the  last  two  expressions  become  equivalent  to  delta  function  impulses  in  the  limit  of 

Ln  ->0. 

Once  the  simplified  mode  shapes  of  equation  29  are  obtained,  the  initial  conditions  for 

v,  w,  q,  and  r  are  used  to  complete  the  solution.  Equations  49  and  50  are  the  analytical  solutions 

for  the  fixed  plane  translation  velocities  v  and  w ,  expressed  in  phase-amplitude  form. 
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(Os-OfXP-Os)Ln 


_DP. 


n{  =  ^P  —  Of  jsn  +  <|>B  +  Of  s  ,  Qs  =-^P-Osjsn-(t>B-<Dss; 
Af  =(p-Of)Ln+Of  ,  As  =  ^P  -  Of  ^Ln  -  Qs . 


(56) 


(57) 

(58) 

(59) 

(60) 


Swerving  motion  is  measured  along  the  earth-fixed  Jj  and  Kj  axes.  To  an  observer  standing 
behind  the  gun  tube,  these  axes  are  oriented  such  that  positive  Jj  is  to  the  right  and  positive  Kj 
is  pointed  downward.  The  swerving  motion  results  from  a  combination  of  the  normal 
aerodynamic  forces,  as  the  projectile  pitches  and  yaws,  plus  the  forces  and  moments  due  to  the 
applied  impulse.  Differentiating  equations  22  and  23  with  respect  to  nondimensional  arc  length, 
generates  the  swerve  equations  such  that 


13 


and 


■\ 

_y 

D 

J 


sgn  s-sn  -  sgn  s  -  sn 


2L„ 


Av 

V0 


(61) 


ft 


> 

z 

=  1} 

sgn(s-sn^ 

!-sgn(s-sn  -Ln) 

)sin(ps  +  <|>B) 

<D  j 

2Ln 

+ 


G-Aw 

V0 


(62) 


Integrating  these  equations  is  straightforward,  but  the  solutions  are  judged  to  be  too  long  to  be 
listed  here.  However,  the  asymptotic  limit  of  these  solutions,  s  ->  oo ,  is  of  special  interest. 

For  a  stable  projectile,  the  swerve  caused  by  epicyclical  vibration  decays  as  the  projectile 
progresses  downrange  and  does  not  affect  the  long-term  lateral  motion  of  the  projectile.  When  a 
lateral  pulse  is  applied  to  the  projectile  at  arc  length  sn ,  its  effect  on  the  target  impact  point  is 

predominantly  due  to  induced  jump,  provided  the  target  distance  is  sufficient  to  allow  the 
transients  to  decay.  Projectile  linear  theory  shows  that  the  long-term  center  of  mass  solution,  or 
swerve,  contains  terms  that  remain  bounded  with  arc  length  s  plus  terms  that  are  linear  with  s , 
and  if  gravity  is  included  the  solution  will  have  even  higher  order  diverging  terms.  These  higher 
order  terms  are  typically  denoted  as  gravity  drop.  The  linear  terms  are  called  jump  terms,  which 
are  caused  by  initial  conditions  at  the  gun  muzzle,  lateral  pulse  forces,  and  aerodynamic 
characteristics.  Mathematically,  setting  gravity  to  zero  and  subsequently  evaluating  the 
following  limits  formally  defines  total  jump. 


lim 

s-*» 


y(s) 

Ds 


lim^7=rK 

s_*oo  S 


(63) 


The  total  jump  vector  T  is  expressed  as  the  sum  of  two  vectors.  The  first  vector  represents  the 
jump  attributable  to  muzzle  conditions,  and  the  lateral  pulse  is  represented  by  the  second,  thus 
resulting  in  the  following  expression: 


_ A 

(af+b)2+(ae+c)2 


Jaf+b 

(AE  +  C 


-AE-Cllv0F  +  w0E-q0Dl  Icosll 
AF  +  B  1 1  w0F-v0E-r0D  I  jsinll 


(64) 


for  which 
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a(af+b)m*+a(be-cf)f* 

(af+b)2+(ae+c)2 

„  (abf+ace+c2+b2)f,-a(ae+c)m* 

(af+b)2+(ae+c)2 

n  =  tan-,(Z/Y)-^s--Psn-cpB  ->  tan^Z/Yj-Ps. -<pB  asLa->0, 


and 

VFTFjgfe#  as  L  ->  0.  (65) 

(PL./2) 

The  quantities  A  and  II  are  the  magnitude  and  phase  angle  of  the  jump  vector  attributed  to  the 
lateral  pulse.  It  is  interesting  to  note  that  the  jump  terms  induced  by  the  lateral  pulse  are  not 
dependent  on  the  lateral  state  of  the  projectile  (v(s),w(s),q(s),  and  r  (s)),  thus  this  particular 

contribution  to  jump  is  not  coupled  to  the  projectile’s  angle  of  attack.  Equation  65  extends  the 
work  by  Guidos  and  Cooper  ( 6 )  who  considered  nonspinning  projectiles  subject  to  a  singular 

delta  function  impulse.  Limiting  ^  Ln,P  0  in  equations  64  and  65  produces  expressions 
that  agree  with  their  previous  predictions. 


6.  Lateral  Pulse  Smearing 


In  order  to  better  understand  the  swerve  response  due  to  a  lateral  pulse,  results  for  a 
representative  spin-stabilized  40-mm  projectile  configuration  are  calculated  and  discussed  in  this 
and  the  next  two  sections.  Nominal  values  for  the  aerodynamic  coefficients,  the  projectile 
physical  parameters,  and  flight  characteristics  are  given  in  the  Appendix.  All  results  presented 
use  these  values  unless  specified  differently.  Diversions  from  these  nominal  values  are  clear 
from  the  context  of  the  particular  chart  under  examination. 


Plots  showing  the  effect  of  smearing  are  given  in  Figures  3  and  4  for  the  applied  force 

to 


F*  -  3.838  x  10  5  and  moment  arms  Xt/D  = 


-0.635 


->  +  0.635 


Negative  values  of  X, 


indicate  that  the  application  point  of  the  pulse  force  is  aft  of  the  mass  center  while  positive 
values  indicate  the  pulse  force  is  forward  of  the  mass  center.  To  illustrate  the  smearing  effects, 
the  roll  position,  <|>N ,  of  the  lateral  impulse  force  is  assumed  to  act  primarily  along  the  nonrolling 

Y-axis.  This  means  the  arc  length,  sN ,  corresponding  to  the  center  of  the  pulse  satisfies  the 
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A  Cos  n  X  10 


Figure  3.  Jump  component-j  due  to  lateral  impulse  vs.  impulse  duration. 


Figure  4.  Jump  component-k  due  to  lateral  impulse  vs.  impulse  duration. 
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expression  «|)N  =  <|)'sN  +  <|>B  =  27iN ,  N  =  0,1,2  •  •  • .  To  ensure  this  force  is  acting  nominally  along 
the  nonrolling  Y-axis,  the  activation  point  begins  at  sn  =  sN  -  Ln  /2  so  that  the  duration  of  the 
impulse  brackets  <|>N  =  2nN .  Increasing  pulse  length,  Ln ,  causes  the  jump  components  to 
cyclically  decay  while  its  value  at  Ln  =  0  corresponds  to  a  lateral  impulse  that  is  proportional  to 
the  delta  function  8(s  -  sn).  Values  of  Ln  where  the  jump  is  zero  represent  situations  where  the 

duration  of  the  lateral  pulse  coincides  with  a  roll  cycle.  This  accounts  for  the  post  multiplier, 

f  \ 


sin 


pDLn 

2Vn 


,  in  equation  65  having  2 % 


zero  crossings  and  causes  the  cyclical  jump  results  in 


Figures  3  and  4.  Notice  for  the  cases  presented  the  response  due  to  a  pulse  in  the  Jj  direction  is 
predominantly  along  Jj ;  however,  a  smaller  out-of-phase  component  of  swerve  is  also 


generated. 


7.  Lateral  Pulse  Response  Magnitude  and  Phase  Angle 


Figure  5  has  charts  showing  the  absolute  value  of  A  as  a  function  of  the  pulse  length  Ln  for  the 

to 


-5 


previous  parameter  values  F  =  3.838  x  10  and  Xr/D  = 


-0.635 


->  +  0.635 


Here  again 


the  impulse  is  assumed  to  bracket  <f»N  =  27iN .  The  maximum  value  occurs  at  Ln  =  0 ,  which  is 

the  delta  function  result,  as  is  justified  from  equation  65.  Consideration  of  A  as  a  function  of 

to 


spin  rate  p  is  displayed  in  Figure  6  for  several  pulse  durations  Ln  = 


0- 


->40 


.  This  figure 


shows  that  the  rate  at  which  the  aerodynamic  jump  magnitude  decreases,  with  spin  rate,  is 
strongly  dependent  on  the  pulse  duration.  Again,  this  is  attributed  to  cancellation  (i.e.,  smearing) 
effects,  which  become  more  pronounced  with  increasing  values  of  Ln . 


Figure  7  gives  the  phase  angle,  II ,  (equation  65)  as  a  function  of  spin  rate  for  the  same  values  of 
the  moment  arm  length  Xr/D  previously  discussed.  Notice  that  the  phase  angle  linearly 

increases  with  spin  rate,  indicating  that  spin-stabilized  projectiles  are  most  susceptible  to  an 
out-of-phase  swerve  response  due  to  lateral  pulse  forces. 


8.  Target  Interception 


The  last  topic  discussed  in  this  report  examines  the  relationship  between  a  given  change  in 

aerodynamic  jump,  J  l  ?  an(j  the  impulse  needed  to  create  this  change.  This  may  be  desirable 
[ATkJ 
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Jump  Magnitude  Ax  10  H  Magnitude  A  x  10 


5.  Magnitude  of  jump  due  to  lateral  impulse  vs.  impulse  duration. 


Figure  6.  Magnitude  of  jump  due  to  lateral  impulse  vs.  spin  rate. 
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Xi/D  =  -0.636 


F  =3.838x10  &  tan  Z/Y-n  =  2jtN 


Xifl>  =  0.0 

/Xr®  =  -0.318 


Xi/D  =  +0.635' 


XitT»  =  +0318 


Spin  Rate  p  (1/s) 


Figure  7.  Phase  angle  of  jump  due  to  lateral  impulse  vs.  spin  rate. 


when  the  activated  impulse  is  strong  enough  to  force  a  projectile  to  strike  its  target,  located  at  a 
relative  position  j^J } ,  with  a  single  impulse.  Hence,  equation  64  requires 


AF,  1.1  sinll  1 
AT,  |  A 1  cosn  |  ‘ 


Solving  this  system  for  II  and  insisting  that  II  remains  real  valued  yields  the  following  two 
equations  for  target  interception: 


•  2  PDs»  pDs„(Arj  +  ArK 

kin  _  —  ,  \ 


2V0  +  L, 


255a. +  *„  =-P^-  +  tan-,f^^— ^^|  +  27tN 


ArT  Sr + Arv  Sr 


Using  equation  65,  the  impulse  force  F,  can  be  determined  from  equation  67,  for  a  given 
moment  arm  Xr ,  and  then  equation  68  solved  for  <|)B  will  guarantee  target  impact,  provided  the 

f  AFj  ] 

relative  target  position,  <  > ,  is  known. 

LArKj 


9.  Conclusions 


An  analytical  approach  for  quantifying  the  effect  of  a  lateral  square  impulse  disturbing  a 
projectile  during  free  flight  has  been  presented.  All  of  the  analysis  was  based  on  projectile  linear 
theory,  which  produces  simple  closed  form  solutions  for  the  assumed  square  pulse  disturbance. 
These  solutions  are  then  used  to  calculate  the  projectile  swerving  motion,  so  that  the  long-term 
effects  of  the  lateral  impulse  are  readily  determined.  The  primary  interest  regarding  target 
interception  is  the  projectile’s  aerodynamic  jump.  Changes  to  aerodynamic  jump  caused  by  the 
lateral  impulse  forces  were  shown  to  produce  easy-to-understand  additive  contributions  to  the 
usual  aerodynamic  jump  of  a  free-flight  projectile  with  no  applied  impulse. 

The  important  question  concerning  smearing  effects  originating  from  a  finite  length  pulse  has 
been  addressed  where  pulse-induced  control  authority  is  shown  to  diminish  as  1/Ln  . 

Calculations  also  show  the  additional  aerodynamic  jump  magnitude,  2 ,  decreases  with  spin  rate 
p  when  subjected  to  various  pulse  lengths  Ln . 
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Appendix  A.  Flight  Coefficients  for  a  40-min  Projectile 


The  numerical  values  used  for  the  graphical  presentations  given  in  this  report  are  shown  in  the 
following  matrices: 

Aerodynamic  coefficients: 


^XO 

0.279 

o 

X 

2.672 

z 

U 

- 

>  «  4 

2.329 

Cypa 

-0.295 

^LP 

-0.042 

^MQ 

-1.800 

Physical  parameters: 


r 

m 

0.0116Slug 

lx 

2.85  xlO'5  Slug  ft.2 

Iy 

2.72  xlO'5  Slug  ft.2 

•D 

►  ~  < 

0.137  ft. 

^LCOp 

0.237ft. 

0.239ft. 

slcg 

0.0713ft. 

(A-l) 


(A-2) 
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Flight  characteristics: 


p  2.38xl0_3  Slug  ft. 

i-»  >  V0  ►  =  <  250.0  ft.  sec.  * 
p  399.7  sec.  ' 
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List  of  Symbols 


C, 

D 

Fa 

F* 

g 

G 

lx 

Iy 

L 

M 

N 

m 

P 

q 

7 

S 

s 

sgn(x) 

u 

V 

w 

V0 

u 

V 

w 


Projectile  aerodynamic  coefficients 

Projectile  characteristic  length  (diameter) 

Dimensional  impulse  force 

Nondimensional  impulse  force  F*  =  DFd/mV^ 

Gravitational  constant 

Scaled  gravitational  constant  G  =  gD/V0 

Mass  moments  of  inertia 

Applied  moments  about  projectile  mass  center  expressed  in  the  no-roll  frame 
Projectile  mass 

Angular  velocity  components  vector  of  projectile  in  the  no-roll  frame 

Surface  area  S  =  7iD2/4 
Nondimensional  arch  length 

jo  T  =  0 

Mass  center  velocity  components  in  the  body  reference  frame 
Forward  velocity  of  projectile 

Mass  center  velocity  components  in  the  no-roll  reference  frame 
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-€  ®-e-  a  >  > 


X. 


p 


Dimensional  moment  arm  length 

Applied  force  components  in  the  no-roll  reference  frame 

Position  vector  of  body  center  of  mass  in  an  inertial  reference  frame 

Longitudinal  aerodynamic  angle  of  attack 
Lateral  aerodynamic  angle  of  attack 

K  Components  of  aerodynamic  jump  due  to  lateral  impulse 
Phase  angle  of  the  aerodynamic  jump  due  to  lateral  impulse 
Euler  roll,  pitch  and  yaw  angles  of  the  projectile 
Euler  roll  angle  of  the  applied  impulse 
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